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Abstract 

The strange quark polarization puzzle, i.e. the contradiction between the 
negative polarized strange quark density obtained from analyses of inclusive DIS 
data and the positive values obtained from combined analyses of inclusive and 
semi-inclusive SIDIS data using de Florian et. al. (DSS) fragmentation functions, 
is discussed. To this end the results of a new combined NLO QCD analysis of the 
polarized inclusive and semi-inclusive DIS data, using the Hirai et. al. (HKNS) 
fragmentation functions, are presented. It is demonstrated that the polarized 
strange quark density is very sensitive to the kaon fragmentation functions, and 
if the set of HKNS fragmentation functions is used, the polarized strange quark 
density obtained from the combined analysis turns out to be negative and well 
consistent with values obtained from the pure DIS analyses. 
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In the absence of neutrino reactions on a polarized target, the inclusive polarized 
deep inelastic lepton-hadron reactions determine only the sum of quark and anti-quark 
polarized parton density functions (PDFs), Aq(x) + Aq(x), and provide no information 
at all about the individual polarized anti-quark densities. All analyses of the polarized 
inclusive (DIS) data have produced results for the polarized strange quark density 
function, As(x) + As(x), which are significantly negative for all values of x (for more 
recent analyses see [H [2])- One way to determine polarized quark and anti-quark 
densities separately is to use the data on polarized semi-inclusive reactions (SIDIS) 
like I + p — > I + h + X, where h is a detected hadron. In the past few years more data 
on polarized SIDIS processes have become available and have led to assertions that 
As(x) + As(x) is positive for most of the range of measured x. In the following we 
discuss possible resolutions to this puzzling state of affairs. 

It should be noted that in the study of the SIDIS data it is usual to simplify the 
analysis by taking As(x) = As(x). On the other hand, it has been suggested, on the 
basis of theoretical models, that this equality is badly broken j3j , and that this could 
be the cause of the conflict. However, it is crucial to realize that: (i) the DIS result for 
As(x) + As(x) is independent of the relationship between As(x) and As(x), and (ii) 
that the COMPASS estimate [I] of the difference As(x) — As(x) is much smaller than 
the theoretical model estimates and thus cannot be the cause of a serious error in the 
extraction of As(x) + As(x) from a combined analysis of the DIS and SIDIS data. 

The key to resolving the puzzle lies, we believe, in the properties of the fragmen- 
tation functions (FFs) needed in the theoretical expressions for the measured SIDIS 
cross-sections and asymmetries, which involve convolutions of either unpolarized or 
polarized PDFs with the FFs. There are three modern versions of the FFs in the 
literature, Hirai et al. (HKNS) [5], de Florian et al. (DSS) [6 J and Albino et. al. 
(AKK) [7], sometimes differing significantly from each other. They are based mainly 
on semi-inclusive e + e~ annihilation data (HKNS), e + e~ annihilation and RHIC data 
on reactions like pp — > tt or K + X (AKK), and a global analysis (DSS) of the data on 
semi-inclusive e + e~ annihilation, the proton-proton collisions at RHIC and unpolarized 
SIDIS processes. 

The early claim by the HERMES Collaboration [8] that the polarized SIDIS data 
implied marginally positive As(x)+As(x) in the measured x range [0.023-0.3] was based 
on a LO analysis of the data. In 2008, de Florian, Sassot, Stratmann and Vogelsang 
(DSSV) carried out a combined NLO QCD analysis [5] of polarized DIS, SIDIS and 
RHIC data using the DSS fragmentation functions and effectively confirmed the LO 
result. More precisely, using the assumption As(x) = As(x) they obtained a sign- 
changing solution for As(x) + As(x), negative for x < 0.03 and positive in the region 
x > 0.03. Later we repeated this analysis [TO] , using polarized DIS and SIDIS data and 
found substantial agreement with DSSV. We confirmed the sign changing behavior of 
As(x), though our As(x) is less negative at x < 0.03 and less positive for large x and 
compatible with zero within the errors. Note that the polarized pp data from RHIC are 
not important for the determination of the polarized quark and anti-quark densities; 



1 



they constrain mainly the gluon polarization. 

After convincing ourselves that the puzzle could not be resolved by taking As(x) ^ 
As(x) we began to try to test whether the problem lay in the properties of the FFs. 
Now the largest disagreements between the various sets of FFs in the literature occur 
in the kaon production sector. To this end we first carried out a combined NLO QCD 
analysis [11] of the polarized world DIS data [12] and just the pion SIDIS data [U [13] , 
using the DSS FFs. Note that in this case only the sum x(As + As)(x,Q 2 ) can be 
determined from the data because of the reasonable assumption = used for all 
the sets of the fragmentation functions. The result for x(As + As)/2 is illustrated in 
Fig. 1 and compared to those obtained from the LSS'06 DIS analysis [2] (red curve) 
and the combined LSSTO fit to the DIS and SIDIS data [10] (black curve). As seen 
from Fig. 1, in the presence of only the A\ N data, x(As(x) + As(x))/2 (blue curve) is 
still negative in the measured x region as in the analyses of the purely inclusive data. 
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Figure 1: Comparison between polarized strange quark densities obtained from differ- 
ent kinds of NLO QCD analyses (see the text). 

This definitely seemed to point towards the kaon FFs as the source of the conflict. 
Note also that it had already been pointed out by COMPASS Collaboration (2nd and 
3rd refs in [T3]) that in the LO QCD approximation the value of the first moment of 
As(x) in the measured range of x is very sensitive to the assumed value of the ratio 
of the s-quark to w-quark fragmentation functions into positive kaons. Therefore, we 
carried out a new combined NLO QCD analysis of the polarized DIS and all the SIDIS 
data [H [T31 E] using the HKNS set of FFs [5], which differ significantly from the DSS 
ones in the kaon sector, especially for the transition s — > K + , as shown in Fig 2 [15]. In 
Fig. 2 two error bands for the HKNS FFs are presented. The narrow one corresponds 
to Ax 2 = 1 while the wide corridor corresponds to Ax 2 = 19.2. The latter value 
corresponds to 17 parameters fit in the MINUIT-procedure when only the statistical 
errors are taken into account. However, the authors of [5] apply this procedure for the 
statistical and systematic errors added in quadrature which definitely overestimates 



1 I — 1 ' ■ — 

LSS'trafDIS) 


XAS=XAS 


" LSS"I0(DIS+SIDIS) 




DIS+SIDIS(onlyA" 1N ) 




' Q 2 = 2.5 GeV 2 

i i i i 


x(As+As)/2 . 



2 



the uncertainties. 
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Figure 2: Comparison between NLO HKNS and DSS kaon FFs at Q 2 = 10 GeV 2 . 



The method used is the same as in our previous analysis [10] of the same set of 
data when the DSS FFs were used. Note that the present SIDIS data are not precise 
enough to determine separately As(x) and As(x). So, as in our previous analysis 
the assumption As(x) = As(x) was used. A good description of the SIDIS data 
(x^ rP =0.92) is achieved using the HKNS FFs (NrP is the number of corresponding 
experimental points). The quality of the fit to the data is demonstrated in Fig. 3 (black 
curves) for some of the SIDIS asymmetries obtained by the HERMES and COMPASS 
Collaborations. The new curves are compared to our previous theoretical curves (red 
ones) obtained from the best fit to the data using the DSS FFs (xjy rP =0.87). As seen 
from Fig. 3 the results from both the fits are very close to each other and for some of 
the asymmetries the curves are almost identical. 

Let us discuss the impact of the HKNS fragmentation functions on the polarized 
sea-quark densities. It is known that the present SIDIS data do not influence the gluon 
polarization. It is mainly determined from inclusive DIS and semi-inclusive pp RHIC 
data. The new values of the sea quark and gluon polarized densities (black curves) are 
presented in Fig. 4 together with their error bands and compared to those obtained 
using the DSS FFs (LSS'10). As seen from Fig. 4 the changes in the polarized sea 
quark densities are as follows: negligible for xAd(x), visible for xAu(x) at x > 0.03 
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Figure 3: Comparison of our NLO LSS'll (black curves) and LSS'10 (red curves) 
results for the SIDIS asymmetries with the data at measured x and Q 2 . 



and dramatic for xAs(x), although the central values of the first moments of As(DSS) 
and As(HKNS) are very close to each other (-0.052 ± 0.016 and -0.048 ± 0.012 at 
Q 2 = 1 GeV 2 for DSS and HKNS FFs, respectively) and coincide within the errors. In 
Fig. 4 our LSS'06 result j2] for x(As(x)+As(x))/2 (blue curve) obtained from the NLO 
QCD analysis of the world inclusive DIS data is presented too. We find now that if 
the HKNS FFs are used, As(x) is negative and well consistent with (As(x) +As(x))/2 
obtained from the pure DIS analyses [U [2] . 

In conclusion, we have found that in the presence of semi-inclusive DIS data the 
strange quark density is very sensitive to the choice of the FFs. We have also demon- 
strated that the strange quark polarization puzzle can be resolved by using the HKNS 
set of fragmentation functions rather than the DSS ones. Finally, we like to stress we 
do not claim to have presented a unique resolution to the strange polarization puzzle. 
Our analysis illustrates only how badly we need to have a more reliable determination 
of FFs in order to extract correctly the polarized sea quark densities. To this end 
precise unpolarized SIDIS cross-section data are very important. 
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Figure 4: Comparison between NLO LSS'11(HKNS FFs) and LSS'10(DSS FFs) sea 
quark and gluon polarized PDFs at Q 2 = 2.5 GeV 2 . The blue curve corresponds to 
x(As(x) + As(x))/2 obtained from the pure DIS analysis [2J. 
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